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ABSTRACT 

O We study the strong coupling limit of the Bethe ansatz solutions in the 

massive Thirring model. We find analytical expressions for the energy eigen- 
values for the vacuum state as well as n— particle n— hole states. This formula 
is compared with the numerical results and is found to achieve a very good 
agreement. 

Also, it is found that the 2— particle 2— hole and higher particle— hole 
states describe n— free bosons states in this limit. The behaviors of the 
strong coupling limit of the boson mass for various model calculations are 
examined. We discuss an ambiguity of the coupling constant normalization 

CN ! due to the current regularization. 
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1 Introduction 

Recent calculations for the massive Thirring model have presented a debate 
over the energy spectrum of the bound states [1-4] . Several different methods 
give different results on the spectrum of the bound state. For a long time, 
people have believed that the semiclassical results by Dashen et al.[5] are 
exact in spite of the fact that they took into account only the lowest order 
quantum fluctuations in the path integral method. However, the recent cal- 
culation by the light cone procedure shows that there is only one bound state, 
and the spectrum of the bound state energy as the function of the coupling 
constant is different from the semiclassical result [1-3]. Further, the recent 
calculations based on the Bethe ansatz solutions [4] present a numerical proof 
that there is only one bound state, and the spectrum seems to be consistent 
with the light cone results. 

In this paper, we present analytical calculations of the strong coupling 
limit of the Bethe ansatz solutions for the massive Thirring model and show 
that the analytical expressions obtained here agree very well with those calcu- 
lated by numerically solving the Periodic Boundary Condition (PBC) equa- 
tions of the Bethe ansatz solutions [6,7]. Here, we obtain the energy eigen- 
values of the vacuum, lp— lh states (symmetric and asymmetric cases ) and 
2p— 2h and higher particle hole states (symmetric case). The analytical for- 
mula shows that the n particle— n hole state is just n times lp— lh state 
energy (a boson mass). Therefore, there is only one boson state in the mas- 
sive Thirring model. This shows that the n particle— n hole states are all 
scattering states. 

Further, we show the behaviors of the strong coupling limit of the bo- 
son mass for various model calculations. It turns out that the analytical 
expression of the boson mass at the strong coupling limit with the Bethe 
ansatz solutions is different from the light cone prediction. This may indi- 
cate that the normalization ambiguity of the coupling constant due to the 
fermion current regularization in the massive Thirring model is different from 
the massless Thirring model. For the massless Thirring model, Klaiber [8] 
proves that the coupling constant has a normalization ambiguity which arises 
from the fermion current regularization. In the case of the massive Thirring 
model, it is expected that the same type of the coupling constant ambiguity 
may well appear. However, we will see later that the normalization ambiguity 
of the coupling constant is more complicated than expected for the massive 
Thirring model. 

In the next section, we briefly describe the Bethe ansatz method which 
is applied to solving the massive Thirring model. In section III, we discuss 
analytical expressions for the energy eigenvalues with the strong coupling 
expansion. Then, section IV treats numerical calculations of the PBC equa- 
tions in the strong coupling region and we compare them with the analytical 
expressions. In section V, we examine the boson mass at the strong coupling 
limit. Finally, section VI summarizes what we have understood in this paper. 



2 Massive Thirring model and Bethe ansatz 
solutions 



The massive Thirring model is a 1 + 1 dimensional field theory with current 
current interactions [4,11], and can be solved by the Bethe ansatz method 
[4,6,7]. In this case, the eigenvalue problem is reduced to solving the following 
PBC equations, 
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where $ denotes the rapidity of the particle and rij's are integer. m , L 
and gs denote the bare mass, the box length and the coupling constant, 
respectively. 

2.1 Vacuum state 

The PBC equations for the vacuum which is filled with negative energy par- 
ticles ( (3i — iir — OLi ) can be written as 
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(i = 0,±l,..,±JV o ), 

where Lq is defined as Lo = m^L and JVo = (N — l)/2. In this case, the 
vacuum energy E 1 ^ can be written as 
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2.2 lp - lh state 

For one particle-one hole (lp — lh) states, we take out one negative energy 
particle (io-th particle) and put it into a positive energy state. In this case, 
the PBC equations become 
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where j3i can be a complex variable as long as it can satisfy eqs.(2.10). These 
PBC equations determine the energy of the one particle-one hole states which 
we denote by E^{ h , 
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2.3 2p - 2h states 

For two particle- two hole (2p — 2h) states, we take out the %\— th and the 
%2— th particles and put them into positive energy states. The PBC equations 
for the two particle-two hole states become 
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In this case, the energy of the 2p — 2/i states E\ J^ becomes 
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Here, we note that the symmetric case ( i\ = —%2 ) always gains the energy 
and therefore is lower than other asymmetric cases of 2p — 2h states. Higher 
particle-hole states are constructed just in the same way as above. 



3 Strong coupling expansion 

Here, we present the strong coupling expansion of the PBC equations. We 
take the limit of gs — ► oo and obtain the energy eigenvalues analytically. 
This limit of gs — > oo can be taken since there is no coupling constant 
renormalization in the massive Thirring model [15]. 

3.1 Vacuum state 

First, we consider the strong coupling limit for the vacuum [16]. We assume 
that gB is much larger than any of the rapidity ctj, namely, 

^B>«i- (3.1) 



In this case, eq. (|Q|) becomes in terms of bi = ^gBL ai, 

Further, the vacuum energy is written as 
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Since frj's have the symmetry of bi = — b_i, we can rewrite eq. (|3.2| ) as 
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From this equation, we can easily obtain 
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Therefore, the vacuum energy can be explicitly written up to 1/gs order, 

E v = -(2N + l)m - — [8N (N - 1) + 12JV ] • (3.6) 
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3.2 lp — lh state ( symmetric ) 



Next, we treat the lp — lh states [15]. We assign the positive energy particle 
by Iq. In the strong coupling limit, eqs.(2.10) become 
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where e(a) denote the step function. These equations have two solutions, 
the symmetric and the asymmetric solutions. For the symmetric case, one 
easily sees, since aij = — a_. 
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Also, for other ctj's, the equations can be rewritten using fej's, 
^ = 16Ett3T, + 4 (i = l,...,iVb). 

In this case, we can evaluate the energy of lp — lh symmetric case in the 
same way as the vacuum and obtain 

E® h = -(27Vo - l)ro„ - -^ [8iV (iV - 1) + 4iV ] . (3.10) 

Therefore, the lp— lh energy for the i$ = case with respect to the vacuum 
becomes 
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3.3 lp — lh states ( asymmetric ) 

Now, we discuss the asymmetric solutions. In this case, we obtain the fol- 
lowing coupled equations up to 1J9b- 
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From the numerical analysis, we can put 

lAol > l«i|- 

Therefore, the above PBC equations are reduced to 
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The energy of the one particle-one hole states E\°[ h becomes 
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Therefore, we should calculate the cosh/? io and the sum of the af. First we 
evaluate the j3i from eq. (|3.13a|) , which can be reduced to a cubic equation 
with x = tanh(/3j /2), 
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It is easy to solve the above equation, and we obtain 
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Next, we calculate the sum of the a?. From eqs.( p.l3b]) and (|3.13c|) , we 
obtain 
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Also, we have from eqs.( p.l3b|) and ( |3.13c| ) when i > i and j < i 
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Therefore, we find 
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Finally, we obtain for the E^°l h , 
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3.4 2j9 — 2h and higher particle— hole states 

Now, we consider 2p — 1h and higher particle— hole states. In this case, the 
symmetric solutions always gain the energy. Therefore, we only treat the 
symmetric solution here. Due to the symmetry, we can easily find 

a = 0. (3.21) 

Let us first consider the 2p — 2h configuration. We assume that the i$— th and 
the — iQ— th particles are in the positive energy state. In this case, eqs.(2.12) 
reduce to 
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Note that this leads to the string— like configurations since the solution for 
b{ becomes pure imaginary. That is, 

b ±l0 = ±2i. (3.23) 

In terms of /3i , this becomes 

P±io = ±"7^- (3-24) 

Therefore, the rapidity interval A of the string becomes 

A = - r t==i + 0\A r \. (3.25) 
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where we explicitly write the behavior of the next order of 1/gs expansion. 
On the other hand, Bergknoff and Thacker [6] assume that the rapidity 
interval of the string behaves for large gs as 

4 
A BT = —i. (3.26) 
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This behavior as the function of gs does not agree with eq. ( |3.25| ) which is a 
solution of the PBC equation. Therefore, the string configurations taken by 
Bergknoff and Thacker are not consistent with the string-like solution that 
satisfies the PBC equations. But this is not at all surprising if one considers 
the way of obtaining the string configurations by Bergknoff and Thacker. 
They assume that the wave functions of the particles should not diverge at 
Xi = — oo as a sufficient condition. However, this cannot happen due to the 
two reasons. The first reason is that one constructs the field theory in the 
box of < Xi < L. Therefore, the boundary is always periodic, that is, the 
wave functions at Xi = L and X{ = are the same. The second reason is 
more physical. The interactions between particles considered here are always 
repulsive. Therefore, the wave functions cannot diverge at any points of the 
space, since they are in the scattering states as bare particles. 

To avoid the confusions, we clarify the string picture. In the nonlinear 
Schrdinger model, the string corresponds to the bound states of the particles 
since they make bound states due to the attractive 5— type interaction. 
However, this is only possible for the bosonic particles. For fermions, there 
is neither two particle bound state nor three or higher particle bound state 
due to the Pauli principle with the 5— type interaction. In the massive 
Thirring model, therefore, we should not consider the string configuration 
which simulates the many particle bound states. 

Now, in the same way as the vacuum case, we obtain the energy for the 
2p — 2/i state 

4?-2fc = -( 2iV o - 3)m - -i- [8N (N - 1) - 4iV ] , (3.27) 
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where we have ignored those terms which vanish when L — > oo with Nq/L 
finite. Therefore, the 2p — 1h energy with respect to the vacuum becomes 
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For higher p — h states, we can evaluate the energy just in the same way as 
the 2p — 2/i state case. For n-particle-n-hole states, the energy with respect 
to the vacuum can be written as 

A<U = 2nm + ^>. (3.29) 

It is important to find that the np — nh state energy is just n times as large 
as that of \p — Ih state energy, that is, 

A<U = nAE[°l lh . (3.30) 

This shows that the n-particle-n-hole states are composed of n free bosons in 
this limit. This result is consistent with the numerical calculations presented 
in ref.[4]. 
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4 Numerical method 



We solve the PBC equations by the Newton method. The type of equation 
we want to solve can be schematically written as 



G(f) = 



(4.1) 



where / = {f\, f 2 , .., f N ) 
is some function. First, 
eq.(4.1) near f as 



are the N variables that should be determined. G 
we denote some initial solution by /q. We expand 



G(fo) + 



= fo + 5x 
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(4.2) 
(4.3) 



We solve this equation for Sx and put them into eq.(4.1). This leads to a 
new set of /, and we consider / as a new f Q and repeat the same procedure 
until we get some convergent results for /. 

This method has a great advantage over the iteration method proposed 
in ref. [4] , namely it gives a good convergence even for the strong coupling 
region. However, the matrix diagonalization can be possible only for a few 
thousand of matrix dimensions if we have to know all of the eigenvalues. It 
is found that the present calculations agree perfectly with those calculated 
in ref. [4]. 

It should also be interesting to check the accuracy of eqs. (|3.6| ), ( |3.10|) 
and (|3.20| ). In Table 1, we show the comparison of the vacuum energies, 
the lp — lh (symmetric) and lp — Ih (asymmetric) energies between the 
analytical expressions and the numerical calculations with the number of 
particles of N = 1601. As can be seen from the table 1, we find quite 
a good agreement between the predictions of the strong coupling expansion 
and the exact numerical calculations. This indicates that the strong coupling 
expansion is indeed a good approximate scheme. 



Table 1 





Numerical 


Analytical 


E v 


-1807.44 


-1809.08 


p(o) 


-1805.20 


-1806.81 


p(i) 


-1767.40 


-1760.44 



The predictions of 1/gs expansion are compared with the computer 
calculations. Table 1 shows the vacuum energy E v , lp—lh energies 
E[_ lh with the symmetric state and E[_ lh with the asymmetric 
state for g B = 245 with N = 1601 and L = 100. 
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5 Bound states at the strong coupling limit 



5.1 Bound state 

Now, we calculate the bound state of the massive Thirring model at the 
strong coupling limit. In ref.[4], Fujita et al. showed that there is one isolated 
boson state and all the other states are continuum states. Therefore, the lp- 
lh continuum energy should start from the free fermion antifermion mass, 
that is twice the physical fermion mass. 
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Thus, the bound state mass M. can be defined as 
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Finally, from eq^^) and eq.( |5.4|) , we obtain the bound state mass as 
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This result can be compared with the prediction of Fujita and Ogura in 
the infinite momentum frame calculation. 
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where go is the Schwinger type coupling constant [8,9]. As can be seen, they 
are different from each other if we assume qb = go- I n ref.[4], it was shown 
that the Bethe ansatz solutions for several cases of the coupling constant 
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are consistent with those of infinite momentum frame calculation by Fujita 
and Ogura [1] with the identification of qb = go- However, it became also 
apparent that the boson mass calculated by the Bethe ansatz solution starts 
to deviate from the light cone result in the strong coupling region. This 
difference may well be related to the normalization ambiguity of the coupling 
constant in the massive Thirring model. 

5.2 Coupling constant ambiguity 

As Klaiber pointed out long time ago, there is an ambiguity of the coupling 
constant in the massless Thirring model. He proves that the coupling con- 
stant is related to each other depending on the regularization. We briefly 
review it in the appendix A. 

The results of the boson mass at the strong coupling limit for the light 
cone and the Bethe ansatz method indicate that the coupling constant am- 
biguity may well be different from eq.(A.8) for the massive Thirring model. 
From the comparison of the numerical calculations and analytical evaluations 
between the light cone and the Bethe ansatz solutions, we can anticipate the 
following relation between g and the Bethe ansatz coupling constant g B 

W2 a(i-ffi- 

7T 2 B + 



9B = g \— r + V R " I > ( 5 - 7 ) 



where B is a free parameter. In this case, the boson mass of the Bethe ansatz 
solutions at the strong coupling becomes identical to the light cone result. 
If we take B ~ 10, then the agreement between the light cone result and 
the Bethe ansatz solutions become very good for whole range of the coupling 
constant. From eq. (|5.7|) , g B becomes for the small value of go 

g B ^ go- (5.8) 

For the large value of go, we obtain 

gB ^ 0o— o-- ( 5 - 9a ) 

This can be derived if we assume the following relation for s and t 

t = ~(^-l)(s-t), (5.9b) 

where X is given as the solution of the following equation, 



7T 7T 



x iJ^- B + 4 -4)x- 2 -4 = . (5.10) 



Up to now, we do not know any physically simple meaning of choosing the 
fermion current regularization which corresponds to eq. (|5.10|) . Further stud- 



ies of the normalization ambiguity of the coupling constant arising from the 
fermion current regularization would be very interesting since we believe 
that it may well be related to some symmetry which is hidden in the massive 
Thirring model. 
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6 Conclusions 

We have presented numerical calculations as well as the analytical expressions 
of the energy eigenvalues of the vacuum and n particle n hole states in the 
strong coupling regions. It is found that the analytical expressions agree very 
well with the numerical values of the vacuum and lp— lh state energies for 
the large values of the coupling constant. 

From the analytical expressions, we find that the 2p— lh and higher parti- 
cle hole states appear as free boson states in the strong coupling limit. This is 
consistent with the recent proof [13] that the S-matrix factorization assumed 
by Zamolodchikov and Zamolodchikov [14] is violated at the quantum level, 
and therefore, the spectrum predicted by the S-matrix factorization is only 
semiclassical. 

We have also obtained the boson mass M. at the strong coupling limit 
analytically. To compare the present result with other calculations, we write 
here the expressions of the boson mass at the strong coupling limit for various 
calculations. 



o 
IT Tfl 

-Mdhn ^ -z , (6.1a) 

2 go 

M FO ^ 2V2tt—, (6.1b) 
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8 m , 

•Mba ^ , 6.1c 

Kg B 

where .Mdhn denotes the result of the semiclassical calculation by Dashen 
et al. , A^fo is the prediction of the light cone calculation by Fujita and 
Ogura, and the present result is denoted by .Mba- If we identify gs = go, 
then .Mba is different from the light cone calculation. 

At the present stage, we believe that the coupling constant normalization 
arising from the fermion current regularization in the massive Thirring model 
is slightly different from the massless Thirring model. It may well be that 
the regularization ambiguity is related to some hidden symmetry which is 
not clearly understood up to now. 



We would like to thank M. Hiramoto and A. Ogura for helpful discussions 
and comments. 
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A Coupling constant ambiguity of the mass- 
less Thirring model 

Here, we briefly review the normalization ambiguity of the coupling constant 

in the massless Thirring model [8,12,17]. 

For the right mover fermion field ipR, we can express it by the massless 

boson fields <f)R and 4>l as 

^ R ~ e is ^-it4>L ( Al ) 

where s and t are free parameters which satisfy the following constraint 

s 2 -t 2 = 4vr. (A.2) 

The s and t can be expressed in terms of the boson coupling constant (3 as 

Now, the fermion current regularization gives another constraint. For exam- 
ple, Schwinger's regularization (go) which makes the fermion current regu- 
larization only in terms of the space coordinate point splitting implies that 

t=£( 8 -t). (A.4) 



In this case, one obtains the following equation, 

p__ i 

4vr ~ 1 + 33. ■ 



(A.5) 



On the other hand, Johnson's regularization (g) which makes the fermion 
current regularization in terms of the space-time coordinate point splitting 
in a symmetric fashion implies that 

t = JL S . (A . 6) 

In this case, one obtains the following equation, 

P 2 2-f 



4tt 2 + * ' 
Therefore, one obtains the relation between go and g as written 



(A.7) 



2fl 
9o = j±j. (Ai 
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